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icon  track  competes  for  the  tame  measurement.  If  two  or  more 


L  SUMMARY 

*  In  the  real-world  multi-target  tracking  problem,  there  exists  the  possi¬ 
bility  for  many  things  to  go  wrong.  Typical  problems  which  arise 
include:  too  few  tracks  are  formed;  too  many  tracks  are  formed  (false 
tracks);  and  inaccurate  position,  course,  and  speed  estimates  are 
reported.  The  above  difficulties  are  often  the  result  of  incorrect  alloca¬ 
tion  of  data  to  individual  tracks.  Algorithms,  while  estimating  the 
motion  of  a  given  target,  inadvertently  mix  in  clutter  and/or  measure¬ 
ments  from  another  target  In  order  for  correct  allocation  of  data  to  a 
given  track  to  be  made,  one  must  have  an  effective  scoring  formula; 
that  is,  some  means  of  determining  how  likely  a  given  assignment  of 
data  is.  To  be  effective,  a  scoring  formula  must  produce  (on  the  aver¬ 
age)  a  better  score  for  correct  assignments  than  for  incorrect  assign¬ 
ments.  Information  useful  in  the  scoring  process  includes  a  priori  intel¬ 
ligence  data  (such  as  initial  target  locations),  models  of  target  motion, 
models  of  the  transmission  channel,  and  expected  moments  of  dutter 
for  the  sensor  pin  setting  being  used.  Basically,  the  score  is  derived 
from  the  residuals  which  come  out  of  the  processing  of  a  batch  of  data 
with  the  extended  Kalman  filter.  This  is  used  to  evaluate  the  likelihood 
of  potential  tracks.  Although  the  “likelihood”  has  an  intuitive  meaning, 
the  term  is  used  here  to  mean  the  probability  density  function  p(x)  of 
the  trackY.  The  expected  cost  of  a  given  assignment  is  derived  with 
the  theory  of  extremals  being  used  to  obtain  the  expected  cost  of  add¬ 
ing  a  clutter  point  in  a  track.  The  resulting  expected  cost  is  then  shown 
to  behave  in  a  quantitative  fashion  and  this  can  be  visualized  from  a 
geometric  viewpoint.  # 

I  I.  INTRODUCTION 

In  the  last  few  years  a  number  of  approaches  to  the  problem  of 
tracking  multiple  targets  in  a  cluttered  environment  have  been  pub¬ 
lished.  Some  aspects  of  the  problem  that  have  been  considered  in 
some  subsets  of  these  publications  include  the  problem  of  false  alarms 
or  missing  measurements,  track  initialization,  multiple  measurement 
Types  (MSI)  and  target  classification. 

Many  of  the  so-called  multi-target  trackers  described  in  the  open 
literature  really  deal  only  with  the  problem  of  one  target  in  clutter. 

This  hypothesis  of  only  one  target  can  greatly  restrict  the  viability  of 
a  multi-target  tracker  to  sort  out  confused  situations. 

Of  the  trackers  that  have  been  proposed,  and  in  some  cases  imple¬ 
mented,  one  can  see  certain  similarities  and  differences  which  allow 
the  trackers  to  be  grouped  into  certain  classes.  The  grouping  and 
certain  applications  of  these  groups  has  led  to  the  following  thoughts. 

Perhaps  the  most  fundamental  aspect  of  a  tracker  is  how  it  handles 
and  interacts  with  the  data.  The  data  is  after  all  our  handle  on  the  real 
world  and  all  the  information  we  have  about  a  specific  tracking  realiza¬ 
tion  is  contained  in  the  data.  A  second  fundamental  aspect  of  the 
multi-target  tracking  in  dutter  problem  is  that  of  alternate  hypotheses. 

It  is  pottible  to  derive  the  probability  distribution  of  tracks  and  clutter 
points  if  one  carefully  specifies  the  a  priori  probability  base,  i.e.,  the 
probabilistic  target  models,  probabilistic  measurement  models,  etc. 

In  very  simple  cases  where  one  can  get  the  optima!  solution,  one  finds 
this  consists  of  all  possible  configurations  of  the  data  into  the  sets.  In 
each  configuration  each  data  set  represents  a  possible  alternate  track 
or  the  set  of  clutter  points.  A  probability  measure  is  assigned  to  each 
possible  total  surveillance  region  picture.  This  globally  optimal 
approach  is  generally  not  reasonable  for  implementation  and  approxi¬ 
mate  or  suboptimal  approaches  must  be  considered. 

Several  approaches  focus  entirely  on  the  construction  of  the  sets  of 
data  (the  construction  of  feasible  tracks).  If  is  this  latter  philosophy 
that  is  being  addressed  in  this  paper. 

Conceptually,  using  maximum  likelihood  techniques,  various  com¬ 
binations  of  data  are  tried  and  then  "scored”  using  log-likelihood 
functions.  The  best  fit  to  the  targei  model  gets  the  lowest  score  and 
this  is  considered  to  be  one  of  the  tracks  in  the  region  if  no  other  low 


targets  compete  for  the  same  measurements,  several  situations  can 
occur.  These  include  the  possibilities  that  the  two  targets  are  lumped 
together,  one  target  is  rejected,  the  targets  get  mixed  with  track  points 
amigned  to  dutter  and  two  “bad”  tracks  reported,  or  the  case  that  all 
points  are  assigned  to  dutter.  It  is  possible,  though  perhaps  not  nor¬ 
mal,  to  find  situations  where  the  choice  of  the  best  nonovcrlapping 
feasible  track  does  not  correspond  to  a  best  surveillance  picture.  This 
is  quite  easy  to  do  if  the  tracks  overlap  and  compete  for  the  same 
measurements. 

Many  trackers  consider  alternate  hypotheses  as  far  as  assigning  mea¬ 
surements  to  a  track.  However,  once  a  decision  has  been  made  that  a 
measurement  belongs  with  another  group  of  measurements,  this  deci¬ 
sion  is  not  re-examined.  Once  the  decision  has  been  made  to  assign  a 
piece  of  data  to  a  “track,”  that  decision  is  final.  This  is  done  because 
the  system  usually  requires  “an  answer.”  Also,  there  is  always  new 
data  coming  into  the  tracker  allowing  new  hypothesis  tests.  In  addi¬ 
tion,  there  is  a  limited  amount  of  computer  response.  One  could  say 
that  the  hypothesis  testing  is  directed  to  make  a  decision  on  the 
proper  surveillance  picture  or  that  the  tracker  is  “decision  directed.” 

In  the  next  section,  we  will  see  how  effective  scoring  algorithms  are 
developed -ones  which  can  handle  the  alternate  scenarios  posed  above. 
Following  this,  in  Section  3,  a  refinement  to  this  scoring  algorithm  is 
proposed  and  it  is  seen  that  the  average  cost  incurred  for 
measurements  to  tracks  can  be  visualized  from  a  geometrical  view¬ 
point.  In  particular,  it  will  be  shown  quantitatively  that  there  extols 
a  unique  number  of  points  (measurements)  in  a  given  track  that  yields 
an  average  minimum  cost.  Incorrectly  assigning  clutter  points  to  this 
track  and  wrongly  assigning  points  to  clutter  will,  on  the  average, 
increase  the  cost  in  a  well-defined  manner. 


2.  SCORING  ALGORITHMS 

In  order  for  a  correct  assignment  of  measurement  data  to  a  given 
track  to  be  made,  we  must  have  an  effective  scoring  formula,  i 
some  means  of  determining  how  likely  a  given  assignment  of  data  to. 

To  be  effective,  a  scoring  formula  must  produce  (on  the  average)  a 
better  score  for  correct  assignments  than  for  incorrect  assignments. 
Information  useful  in  the  scoring  process  includes  a  priori  intelligence 
data  (such  as  initial  target  locations),  models  of  target  motion,  and 
expected  amounts  of  clutter  for  the  sensor  gain  setting  being  used. 
Basically,  the  score  is  derived  from  the  residuals  which  come  out  of 
the  processing  of  a  batch  of  data  with  the  extended  Kalman  filter.  Tim 
is  used  to  evaluate  the  likelihood  of  potential  tracks.  Although  "Tikeli- 
hood”  has  a  useful  intuitive  meaning,  we  use  the  term  to  mean  the 
probability  density  function  p(  X)  of  the  track  X.  The  concepts  we  use 
are  well-known,  since  most  of  the  work  in  estimation  theory  pertains 
to  situations  where  all  the  observations  Z  *  [zj.  z^,  •  •  •,  zn]  are  due  to 
a  single  target.  An  obvious  example  is  the  stochastic  linear  system 

*k-H  *  Ak*k  +  Vie  k  *  0.  !.•••.».  (I) 

*k  -ck*k**k*  k-l.  •••,».  (2) 

with  states  (xk)  C  R*.  observations  { zk]  c  Hl,  process  noise  {  ujJ  c  RM. 
and  measurement  noise  { vk)  C  Rv.  Ak,  B^,  are  matrices  of  appropri¬ 
ate  dimension  that  may  vary  with  time.  The  initial  state  xq  to  a  Gaussian 
random  vector  with  covariance  Pq,  independent  of  the  processes  (uk) 
and  ( vk ) ,  which  are  themselves  zero  mean  white  Cautsian  none  with 
covariances  {Q|J  am!  [  Rk) .  Under  these  assumption*,  the  well 
known  Kalman  equations  provide  minimum  variance  unbiased  estimates 
{  £k]  of  the  states  based  on  all  past  data: 


Ak*k  *  *t<*k4|'t'k»lAkSk*> 

(i) 

\<WWkckvl*»W''- 

(4) 

VkAlc*«Wk- 

(5) 

(l-KkCk+))?k.  k-0. 

(6) 

When  nonlinear  measurements  arc  involved,  a  simple  linearization 
process  (the  extended  Kalman  filter)  is  used.  The  equations  remain 
exactly  the  same,  except  that  the  term 

zk+l  *ck*lAksk  W 


in  the  first  equation  above  is  replaced  by 

*k+l  'hk*l<Aksk>*  <8> 


where  h^+jC*) denotes  the  nonlinear  relabionship  between  the  measure¬ 
ments  and  state  vector,  and  the  C^-matrix  becomes 


Ck 


(9) 


It  is  natural  to  compute  the  likelihood  function  p(X  )  for  the  track 
AC  Z  based  on  the  Kalman  filter  state  estimates.  The  innovations 
sequence  is  an  integral  part  of  this  computation,  which  is  a  sequence  of 
the  measurement  residuals: 


*k+l  *  *k+l  •ck+lAk*k-  k-0,  (10) 


The  (negative)  log  likelihood  function  it  given  in  terms  of  fi  |J  by 

n  n 

c(X)  -  ndim(z)  In  2»  ♦  1  £  In  I  Vk  |  ♦  A-  23  «ltvk1*k-  <»•> 
"  k-1  "  k-1 


The  covariance  matrix  for  the  measurement  residual,  V^,  can  be  com¬ 
puted  directly: 

Vk  -  WkCk,  ♦  Rk+1.  k  “  1.  2,  n  .  (12) 


Each  feasible  track  is  the  result  of  a  hypothesis  test  that  uses  the 
track  likelihood  function  p(  A)  (or  equivalently,  the  negative  log  likeli¬ 
hood  c(X))  determined  from  the  Kalman  filter.  Since  the  density 
function  of  the  alternative  hypothesis  (that  X  is  not  a  track)  is  unknown, 
the  decision  rule  is  simply 


-c(A)  *  In  p(A  |  fSk  }J|m|)  >  on  -  A  f  F,  (13) 

-c(A)  »  In  p(X  |  {Sk)£.,)  <  on  -X  4  F  .  (14) 


Based  on  the  log  likelihood  decision  function,  the  feasible  track  set  is 
F  -  |x|X  Zmlnp(X|  fSk)J.,)  >  o„|.  (IS) 

Primarily,  the  only  random  component  of  c(  X)  is 

li'v's  ,  (16) 

which  for  real  tracks  is  a  chi-squared  random  variable  with  n  *  dimension 
(z)  degrees  of  freedom.*  Therefore,  erro*  probabilities  can  easily  be  com¬ 
puted  for  the  hypothesis  lest  to  predict  the  accuracy  of  feasible  track  con¬ 
struction.  This  has  a  critical  impact  on  the  ultimate  accuracy  of  the  track¬ 
ing  algorithm,  since  a  real  track  mistakenly  excluded  from  F  Cannot  be 
used  in  the  subsequent  Bayesian  decision  process. 


*Om  io  the  fandom  natwe  of  the  mcaiurvnwnt  atrivat  time  and  of  the  wncor 
wMth  kiwi  "Iht  neat  mcMorcfnt  nt."  the  dvrrtopnwnt  of  Vk  to  a  too  raodOM 
in  Mioft  NHn  Hard  to  compare  in  meaning  from  one  realisation  to  the  wit. 


lie  line  any  new  |*uinl  is  added  to  a  pailial  Hack,  it  vlwHild  |m««  j 
coarse  test. 

IliklL  <  <n» 

9  Tine  t«i. 

qn1  «k  <  ’k  «•» 

and,  finally,  the  likelihood  test.  The  coarse  lest  checks  the  magnitude  of 
the  maximum  component  of  the  vector  6^  f  R2  against  and  is 
included  because  it  is  computationally  cheaper  to  perform  than  the  fine 
test.  The  constants  0^  and  can  be  chosen  so  that 

{zkl4vi'«k<n({*kMI«kH-  <  «k}-  «•**! 

One  difficulty  with  this  approach  is  that  because  of  the  “deterministic” 
terms  in  the  likelihood  function 

n 

n  dm(z)  in  2#  +-L  £  ,n  IVkl  (20) 

*  k*l 

it  is  difficult  to  compare  some  of  the  tracks  of  different  length.  It  is 
also  difficult  to  assess  the  absolute  goodness  of  a  score.  Therefore,  an 
alternate  score  with  a  more  absolute  meaning  can  be  defined.  This  is 
described  in  the  next  section. 

3.  REFINED  SCORING 

For  track  i  at  stage  k  define  the  stagewisc  chi-squared  score 

4"<4*4>Tvkl<4'|i'>  •  <21  > 

This  has  the  fcaturcs-if  all  measurements  have  been  assigned  to  the 
correct  targets  and  all  filter  parameters  chosen  cortectly-that 


ElSL)  -  E{Sk) 

m  2 

122) 

E{(Sk-2)2)  -  4 

(23) 

Q 

■ 

K> 

(24) 

f(S)  ■  e'S/2U(S); 

U(S)»0,  s  <  0 

U(S)«  1,  S  >  0 

(25) 

where  E  (• }  indicates  the  expected  value  operator,  a  is  the  standard 
deviation  of  S^,  and  f(S)  is  the  appropriate  density  function  for  a  two- 
dimensional  random  variable  S .  Define  the  cumulative  chi-squared  score 
as 

ni 

s'  -  £  .  •  CO) 

k*  I  _  _  _______ __ _ 

For  easy  evaluation  on  a  single  track  an  evaluation  cost  that  would  be 
quite  meaningful  would  be 


r  •  — L  s' 

‘  Nj 

(27) 

For  this  cost  function  the  statistical  parameters  are 

E(Ci>  -  2 

(28) 

E  ((Cj  •  2)2)  -  4/Nj  *  2/NiW 

(29) 

t(<y)  -  n  y  1  *'y,~  Wy)  • 

2T(Nj) 

(30) 

For  display  purposes  the  use  Of  Cj  as  a  value  measure  of  a  single  track 
makes  a  great  deal  of  intuitive  sense. 

Values  of  Cj «  2  for  reasonable  length  tracks  tend  to  imply  that  the 
filter  parameters  are  set  too  loose.  Thus,  by  reducing  0*  aml/or  R^,  one 
could  obtain  tighter  tracks.  Tracks  for  which  Cj  »  2  clearly  represent 
bad  data  assignment  which  should  not  he  kept.  More  precisely.  If 


2 


Cj  >  2  ♦  «2/N,W  (Jl) 

when  o*3  (for  a  3-tigma  caw)  indicates  either  that  incorrect  data  has 
been  assigned  to  the  track  or  that  tl»e  filter  parameters  (Q.R)  arc  too 
small.  While 

Cj  <  2  -  o2/Nft,  Nj  >  3  for  a  -  3  (32) 


indicates  that  the  filter  parameters  (Q.R)  are  set  too  large. 

For  a  total  surveillance  region  picture  of  M  measurement  points  and 
L  tracks,  one  has: 


Track  No. 

Score 

Number  of  Points  ! 

m  Track  1 

1 

C| 

", 

2 

C: 

Nj 

L 

cL 

..  _  - - 1 

"L 

The  score  for  the  total  area  should  be  made  up  of  these  scores  and  the 
cost  for  assigning  a  point  to  clutter.  The  number  of  points  in  clutter  is 
N^.  the  number  of  points  assigned  to  targets  is  Np,  and  the  total  number 
of  measurements  are  M.  These  arc  related  by 

i 

M  -  Nc  +  Np  Np  “  E  Ni  • 

j*| 

(33) 

A  meaningful  score  could  be  defined  as: 

L 

s  -  E  Nj  Cj  ♦  Nc  Sc 
i-1 

(34) 

where  S£  is  a  score  defined  for  a  clutter  point.  If  all  the  measurements 
are  correctly  assigned  to  the  track: 

i 

E(S)  “  2  £  N;  +  NCSC 
i*l 

(35) 

«  2(M  -  Nc)  ♦  NcSc  . 

(36). 

If  we  define  the  clutter  score.  Sq  and  a  total  surveillance  score  as 

c  *  JL  s 

M 

(37) 

where,  if  all  the  measurement  points  belong  to  the  track, 

E(C)  -  ^  E<s>  -  T7  (2(M-  Nc)  *  N CSC) 

M  M 

(38) 

,,0.2^  *^5 r 
'  M  M  C 

(39) 

1 


truck  will  project  to  u  pi ven  |Nniii  hi  iiieii%utcMn*«il  quee.  I'«h  esantpK  , 
consider  the  two-dimensional  measurement  vector (r.si described  in 
Figure  2. 

If  a  clutter  point  is  added  lo  the  track,  the  increase  in  wore  on  the 
track  will  be  given  by  Equations  (21 )  and  (37).  Tlie  difference  is  lha* 
now  the  value  of  the  score  can  be  written  as 

(  L  ) 

C  "  M  “  l!r  NjCj*(Nc.|)Sc-Cl  j  .  «■> 


umiuw  all  track 
point i  art  properly 
asMpncJ 

Here  Cf  is  the  cost  for  assigning  one  clutter  point  to  the  track.  C1) 
is  a  random  variable  and  its  distribution  depends  on  the  distribution 
of  the  clutter  points.  From  Equation  (20, 

Cl  -  <4-*k>Vkl<lk-*ik>  <«2> 

where  z£  are  the  r,  a  points  for  the  clutter  point.  The  distribution  of 
this  random  variable  depends  on  the  random  nature  of  the  clutter. 

If  the  distance  of  the  dutter  point  from  the  predicted  point  is  assumed 
to  be  Gaussian  in  r  and  *  with  zero  mean,  the  distance  from  the  pre¬ 
dicted  point  will  be  Rayleigh  distributed  and  the  cost  or  score  ,  the 
weighted  square  of  the  distance  (42),  wilt  have  an  exponential  density 


fCj(c) 


u(c)  *  0, c<0 
u(c)  *  l.  c  >0 


(43) 


where  o“  is  a  measure  of  the  dispersion  of  the  dutter  points  with 
respect  to  the  mid-  or  tracker-predicted  point.  The  uncertainty,  0, 
takes  into  account  the  unequal  variance  in  r  and  o  and  is  given  by: 


(range  of  o  in  surface)  (range  of  r  in  surface) 
p  *  ■  -  ^ 


(44) 


The  score  for  the  closest  clutter  point  to  a  point  predicted  by  the 
tracker  (closest  in  the  sense  of  having  the  smallest  score  defined  by 
(42))  will  be  distributed  as  follows  (x  *  score  of  the  closest  point): 

Fx(x)  *  probability  that  the  closest  clutter  point  out  of 
having  a  wore  less  than  or  equal  to  x. 

Fx(x)  -  1  -  e**/2o*  _  (45) 

From  the  theory  of  extremals ,  it  is  obvious  that  this  distribu¬ 
tion  is  identical  to  the  probability  distribution  of  the  event  defined 
below. 

•  ,  V 

Fx(x)  *  (probability  that  all  of  the  clutter  points  have 
scores  less  than  or  equal  to  x} 

or 

Fx(x)  *  1  *  (probability  that  all  clutter  points  have 
scores  greater  than  x) 


so  that  a  plot  of  the  score  for  a  global  surveillance  region  for  a  real  case 
of  Np  measurements  assigned  to  N*  targets  with  all  points  assigned  cor¬ 
rectly  and  Nc  dutter  points  can  be  geometrically  described  by  Figure  I. 

The  curve  from  N*0,  that  is,  all  points  assigned  to  clutter,  to  N*Np, 
the  correct  number  of  points  assigned  to  clutter,  is  just  a  straight  line 
described  by: 

where 

varies  from  M  to  M-Np . 

Beyond  the  point  (N*Np).  a  dutter  point  must  be  assigned  to  the  track 
(avHimiftii  just  one  track  in  cluticr)-  The  cmf  to  assign  a  single  clutter 
point  lo  the  track  can  be  calculated  in  the  following  manner:  A  given 


From  (43).  it  follows  that  the  probability  that  the  score  of  any  one  of 
the  N'c  clutter  points  being  greater  than  x  is  given  by: 

e**/2o2  .  (44) 

Assuming  the  dutter  points  are  independent  of  one  another,  the 
probability  that  all  of  the  Nq  clutter  points  have  scores  greater  than 
x  is  given  by 

(-Nc*/2o2 


and  the  probability  that  all  of  these  clutter  points  have  scores  less  than 
or  equal  to  x  is: 


l-c 


Ncx/2#3 


(47) 


3 


mm!  density  luwltmi: 


2r 


<4») 


T1»e  expected  value  of  the  cost  of  assigning  the  closest  clutter  point  to  a 
track  is  then  gwen  by: 


E <»)-/  xf(*)  dx  -  Isl.  . 

i.  NC 


<«9) 


The  expected  cost  for  assigning  the  next  nearest  point  to  a  track  is 
*  ^ 


(Nc-I) 


(50) 


and.  finally,  the  expected  cost  for  adding  the  k*h  closest  clutter  point 
out  of  to  a  track  is  given  by: 


lZ2 


N  r+  1-k 


(51) 


Note  that  in  computing  the  expected  value  of  the  total  cost  function 
defined  in  (41).  one  must  also  account  for  the  linear  decrease  in  cost 
caused  by  the  decreasing  weighting  coefficient  (N^-k)  on  Sq  when  k 
clutter  points  are  incorrectly  assigned  to  tracks.  This  amounts  to  a 
decrease  of  S^/M  for  every  clutter  point  assigned  to  a  track.  Therefore, 
the  net  increase  in  the  expected  cost  by  assigning  the  k  closest  clutter 
points  out  of  Nf  to  tracks  is  given  by  : 


In  elicit,  Him  b  llie  sum  ul  an  increasing  liypet Indie  Iuiwimmi  and 
a  decreasing  linear  function. 

Since,  on  the  average,  o*  »&c*  the  shape  of  the  right-hand  side  of 
the  curve  in  Figure  I  is  concave  and  monotonically  increasing. 

Note  the  two  regions  of  this  figure.  In  one  part  fewer  points  are 
assisted  to  tracks  than  actually  are  available,  i.e.(  too  many  detections 
were  missed.  In  the  other  half  of  the  figure,  false  alarms  or  clutter 
points  were  assigned  to  tracks.  In  both  pans  of  the  curve  the  assump¬ 
tion  is  made  that  the  data  association  is  done  in  an  optimal  manner  on 
the  average.  Thus,  if  a  clutter  point  is  added  to  a  track,  it  is  the  clutter 
point  Iha!  lies  “closest**  to  a  track  of  those  urussigned.  The  curve  aho 
assumes  that  all  measurements  correctly  assigned  as  track  points  arv 
assigned  to  the  correct  track.  Wrong  assignments  of  data  would,  of 
course,  make  for  even  worse  scores  “on  the  average.*' 

4.  CONCLUSIONS 

In  the  last  section,  it  was  shown  that  (lie  refined  scoring  algorithm 
possesses  appealing  properties  from  a  geometric  viewpoint.  There 
exists  a  unique  number  of  points  in  a  track  that  results  in  lower  average 
cost  than  any  other  number.  Also,  the  sensitivity  of  the  score  to  varia¬ 
tions  in  assigned  number  of  track  points  can  be  controlled  by  the  clutter 
score.  S£.  This  is  readily  apparent  from  Figure  1. 

This  scoring  algorithm  is  therefore  a  very  useful  approach  in  extract¬ 
ing  clutter  points  out  of  a  givrn  target  track. 

The  algorithm  is  currently  being  applied  to  an  ocean  surveillance 
problem  and  the  results  of  this  are  very  encouraging.  For  a  given  data 
set  having  a  false  alarm  rate  of  1 0"*,  i.e.,  one  clutter  point  in  every  104 
measurements,  using  the  refined  scoring  algorithm  defined  in  Section  3, 
we  have  found  that  we  can  effectively  decrease  this  false  alarm  rate  to 
10-7.  This  represents  a  three-order-of-magnitude  decrease  and,  hence, 
the  detection  capabilities  of  the  tracking  algorithm  have  been  signifi¬ 
cantly  enhanced. 


N  -  NUMSt  A  Of  POINTS  ASSIGNED  TO  TRACKS 


Figure  /.  Geometric wt  description  of  refined  storing  elgorlthm. 
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Figure  2.  Two-dimensional  measurement  space* 


